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摘 要 : 针对 高 光谱 图 像 解 混 问 题 进行 研究 ， 发 现 传统 解 混 算法 在 保持 端 元 数目 不 变 的 情况 下 ， 
(ANN) 提出 一 种 估计 单 像素 点 中 端 元 数目 和 类 别 的 解 混 算法 。 首 先 利 用 人 工 神 经 网 


高 。 为 此 ， 基 于 人 
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得 到 的 解 混 精 度 不 


络 对 站 感 图 像 中 各 个 像素 的 端 元 数目 和 关 别 进行 估计 ， 之 后 依据 估计 结果 确定 解 混 算 法 的 目标 函数 ， 并 引入 改进 的 
差分 搜索 算法 对 目标 函数 进行 优化 求解 ， 最 终 获 取 地 物 丰 度 和 待 求 参数 ， 实 现 高 光谱 图 像 的 解 混 。 仿 真 数据 和 真实 
竺 合 实 际 场景 的 情况 。 


的 遂 感 数据 实验 表明 ， 与 现 有 的 解 混 算 法 相 比 ， 所 提 解 混 算法 具有 更 高 的 解 混 性 能 ， 更 加 符 
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Unmixing of hyperspectral images based on endmember estimation of artificial neural network 


Zhang Heng!, Jia Zhicheng!, Chen Lei? ?, Guo Yanju! 
(1. School of Electronic Information Engineering, Hebei University of Technology, Tianjin 300401, China; 2. School of 
Precision Instrument & Opto-Electronics Engineering, Tianjin University, Tianjin 300072, China; 3. School of Information 


Engineering, Tianjin University of Commerce, Tianjin 300134, China) 


Abstract: Aimed at the problems of hyperspectral unmixing, it was found that the unmixing accuracy of the traditional 


unmixing algorithm was not high when the number of endmember was kept constant in unmixing. Thus, based on the 


artificial neural network (ANN) , this paper proposed a novel unmixing algorithm of estimating the number and category of 


endmember in a single pixel. Firstly, the unmixing algorithm used the artificial neural network to estimate the number and 


category of each mixed pixel's endmember in the remote sensing image. Then, the objective function of the algorithm was 


determined based on the estimation results, and the improved differential search algorithm was introduced to solve the 


objective function. Finally, the abundances and the parameters are obtained to realize the unmixing of hyperspectral images. 


The experimental results on simulated and real hyperspectral data demonstrate that compared with the existing unmixing 


algorithms, the proposed unmixing algorithm has higher performance and is more in line with the actual scene. 


Key words: hyperspectral images unmixing; ANN; endmember estimation; differential search algorithm 


0 ”引言 就 被 遥感 仪器 接收 成 像 ， 一 般 只 适合 线性 混合 的 地 物 场 景 ， 
rt ifj NLMM 则 考虑 了 地 物 之 间 在 传递 过 程 上 的 多 层 散 射 作 用 ， 
高 光谱 图 像 解 混 是 指 从 观测 的 高 光谱 图 像 数 据 中 分 解 得 可 用 于 描述 复杂 的 多 线性 混合 的 地 物 场 景 ,更 具有 实际 意义 。 
到 地 物 的 基本 组 成 成 分 〈 称 为 端 元 ) 以 及 每 种 成 分 在 各 个 像 但 是 现 有 的 基于 像 元 混合 模型 的 算法 无 论 是 线性 的 还 是 
元 中 所 占 的 比例 ( 称 为 丰 度 )。 由 于 受到 空间 分 辩 率 以 及 地 物 非 线性 的 , 通常 都 假设 在 解 混 的 过 程 中 , 端 元 数目 是 不 变 的 ， 
复杂 性 的 限制 ， 高 光谱 遥感 器 获取 的 图 像 通 常 是 由 混合 像 元 即 各 个 像素 点 的 端 元 数目 都 是 相同 的 。 而 实际 遥感 图 像 在 形 
中 组 成 的 。 因 此 ， 如 何 从 混合 像 元 中 提取 端 元 和 丰 度 成 为 国 成 的 过 程 中 ， 图 像 中 的 大 多 数 混合 像 元 只 是 由 整 幅 遥感 图 像 
内 外 众多 学 者 研究 的 核心 问题 站 。 近 年 来 ， 一 些 学 者 通过 建 端 元 的 一 种 或 几 种 组 成 ， 即 某 些 混合 像 元 并 非 所 有 的 端 元 都 
立 数学 模型 来 对 混合 像 元 的 形成 过 程 进行 描述 ， 将 高 光谱 解 参与 了 混合 。 因 此 ， 如 果 在 解 混 的 过 程 中 ， 只 选择 相关 的 端 
混 问 题 转换 为 对 像 元 混合 模型 参数 的 优化 问题 ， 从 而 简化 ] 元 进行 解 混 ， 则 会 相应 的 提高 解 混 的 精度 。 本 文 ， 首 先 对 构 
遥感 图 像 的 解 混 过 程 ， 提 高 了 解 混 的 精度 和 效率 ， 取 得 了 一 建 的 人 工 神 经 网 络 Cartificial neural network, ANNO 加 进行 
定 的 研究 成 果 B5。 训练 ， 之 后 利用 训练 好 的 网 络 在 解 混 之 前 对 每 个 像 元 的 端 元 
于 对 像 元 混合 过 程 抽象 程度 的 不 同 ， 和 常见 的 混合 模型 数目 和 类 别 进行 估计 ， 同 时 根据 估计 结果 以 及 地 物 丰 度 约 束 
有 线性 混合 模型 (linear mixing model, LMM) [SI 以 及 非 线 性 项 确定 各 个 混合 像 元 解 混 的 目标 函数 ， 最 后 采用 改进 的 差分 
混合 模型 (nonlinear mixing model, NLMM) Ul, $, LMM 搜索 算法 〈differential search algorithm, DSA) 四 对 目标 函数 
假设 一 个 光子 只 能 代表 一 种 地 物 ， 且 地 物 之 间 没 有 相互 作用 进行 优化 求解 ， 最 终 实 现 解 混 。 
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1 ”高 光谱 图 像 混 合 模型 
1.1 线性 混合 模型 
LMM 假设 遥感 图 像 中 每 个 像素 的 观测 值 等 于 各 端 元 的 


光谱 特征 按照 对 应 的 丰 度 进行 的 线性 混合 。 其 模型 为 
y-Ma-cn (1) 

Ks ?为 含有 工 个 波段 和 N 个 像素 的 观测 矩阵 ，M 是 光谱 

H SABER. a 为 丰 度 矩阵 ，7 为 白 噪 声 。 其 中 ， 由 于 丰 度 矢量 


a 代表 像素 点 中 各 端 元 所 占 的 比例 ， 其 必须 满足 丰 度 非 负 


( Abundance Nonnegative Constraint ，ANC) 和 丰 度 和 为 一 
( Abundance Sum-to-one Constraint，ASC) 的 约束 条 件 00， 即 


a, >05, =1 2) 

1.2 非 线 性 混合 模型 
双 线 性 混合 模型 U1 (bilinear mixture model, BMM) 是 
近年 来 提出 的 可 用 于 微观 尺度 上 地 物 精确 光谱 分 析 的 一 类 


Tr 
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2 人工 神经 网 络 

人 工 神经 网 络 是 一 种 模仿 人 脑 神经 元 的 行为 特征 ， 进 
分 布 式 并 行 信息 处 理 的 数学 模型 ， 通 过 合理 的 改变 系统 中 
量 节点 〈 神 经 元 ) 之 间 相 互 连 接 的 权重 ， 从 而 达到 对 信息 
行 预测 处 理 的 目的 ， 同 时 具有 较 高 的 容错 性 和 自 适 应 学 习 
力 。 近 年 来 ， 人 工 神经 网 络 的 研究 工作 不 断 深 入 ， 已 经 取 
了 很 大 的 进展 ， 其 在 模式 识别 、 医 学 、 农 业 、 图 像 等 领域 
成 功 地 解决 了 许多 现代 计算 机 难以 解决 的 实际 问题 ， 表 现 
了 良好 的 智能 特性 106-135]。 

人 工 神 经 网 络 有 多 种 模型 结构 ， 典 型 的 有 多 层 感知 机 
(multilayer perceptron, MLP) D% 模型 、Hopfield 模型 201、 
BAM 双向 联想 记忆 神经 网 络 20 等 ， 其 中 ， 由 于 经 典 的 多 层 
感知 机 模型 结构 简单 、 计 算 量 小 、 具 有 高 度 的 非 线性 分 类 性 
能 ， 本 文 基于 MLP 的 算法 理论 ， 搭 建 人 工 神经 网 络 ， 其 结 
构 如 图 1 所 示 。 整 个 网 络 由 输入 层 、 隐 含 层 、 输 出 层 构 成 ， 


S 
i 


[lu gà RB ORE D 


上 上 


线性 模型 ， 其 更 加 符合 高 光谱 图 像 混合 的 实际 形成 过 程 。 


并 且 用 大 量 的 节点 将 其 连接 起 来 ， 图 中 用 圆圈 代替 节点 《〈 神 


型 的 BMM 主要 包括 Fan 模型 03 (fan model, FM), GBM 模 
7803] (generalized bilinear model, GBM), PPNMM 模型 14 
(polynomial post-nonlinear mixing model, PPNMM) 等 。 其 中 ， 
GBM 是 一 个 广义 通用 模型 ， 该 模型 假设 遥感 仪器 接收 到 的 
是 经 过 多 种 地 物 依次 反 向 散射 的 混合 光谱 信号 ,其 定义 如 下 : 

y- Mas Y Y yaam, Om, «n (3) 


其 中 ;yj 为 非 线 性 系数 ,控制 光谱 之 间 的 非 线 性 效应 的 幅 值 。 

当 0<yy<1 时 , GBM 可 以 较 好 的 诠释 光谱 信号 的 非 线性 作用 ， 
当 yy=1 时 ，GBM 就 演变 成 Fan 模型 ， 当 yy=0 时 ，GBM 就 
退化 为 LMM 而 PPNMM 是 2012 年 由 Altmann 等 人 提出 的 
一 种 双 线 性 模型 ， 其 采用 多 项 式 级 数 来 描述 端 元 之 间 所 提供 
的 非 线 性 作用 ， 且 考虑 了 对 像 元 光谱 具有 重要 贡献 的 端 元 自 
身 的 二 次 散射 ， 其 模型 表达 式 为 


y-Ma-«b(Moa)O(Ma)-«n (4) 


Rp: MIIE, a 为 丰 度 矩阵 。 
于 BMM 忽略 端 元 之 间 的 三 次 及 以 上 的 高 阶 散 射 ， 在 

些 地 形 复杂 的 场景 解 混 效 果 不 太 理 想 , 因此 ，Heylen 等 人 
在 2015 年 提出 了 无 穷 阶 非 线性 混合 模型 (multilinear mixing 
model, MLM) UI, HI 


liri 


y=- p)Ma+ 


a- PP% Xaa; (n, Om;)+ 
i=l j=l (5) 


R R R 
d- p)p? Y Yaaa O x; Ox)t- 
i=l j=l k=l 


AL 


中 ， 高 阶 散 射 系数 p 使 得 MLM 在 BMM 的 基础 上 将 混合 
模型 扩展 到 无 穷 阶 的 非 线 性 混合 模型 ，MLM 假设 当 光 线 照 
射 到 某 一 地 物 时 ， 会 以 概率 p 与 它 地 物 进 行 非 线性 作用 ， 而 
以 1-p 的 概率 被 遥感 仪器 接收 ， 同 时 ， 当 概率 p=0 时 , MLM 
退化 成 LMM。 
基于 像 元 混合 模型 的 解 混 算法 通常 是 在 端 元 光谱 已 知 的 
础 上 ， 根 据 丰 度 约束 以 及 误差 最 小 化 原则 确立 目标 函数 ， 

而 将 解 混 问题 转换 为 对 目标 函数 的 最 优化 问题 。 因 此 ， 本 
分 别 基 于 LMM、GBM、PPNMM、MLM, 利用 ANN X13 
图 像 各 个 像 元 参与 模型 混合 的 端 元 数目 和 类 别 进行 估计 ， 

根据 估计 结果 确定 解 混 算法 的 目标 函数 ， 并 采用 收敛 效果 更 


经 元 ), 每 一 个 节点 都 代表 一 种 激励 函数 , 每 两 个 节点 之 间 都 
有 一 个 通过 该 信号 的 权 值 ， 被 称 为 权重 。 通 过 已 知 的 输入 、 
输出 ， 不 断 的 反 向 传播 神经 元 的 误差 (Error) 来 更 新 前 馈 网 
络 的 权 值 参数 和 矩阵， 最 终 获取 最 优 网 络 。 
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图 1 人 工 神经 网 络 结构 示意 图 
Fig. 1 Structural representation of artificial neural network 
3 解 混 算 法 描述 
3.1 ANN 端 元 估计 
本 文 基于 ANN 对 遥感 图 像 各 像素 点 的 端 元 数目 和 类 别 
进行 估计 ， 主 要 包括 三 个 过 程 。 
1) 设置 网 络 的 拓扑 结构 
本 文采 用 的 ANN 输入 层 的 节点 数 为 高 光谱 图 像 的 波段 
数 ， 隐 含 层 的 节点 数 n 为 


n=NN+M +m (6) 


其 中 : N 和 M 分 别 为 输入 层 和 输出 层 的 神经 元 个 数 ; m 为 1 
到 8 之 间 的 整数 ， 输 出 层 节 点 的 数目 等 于 所 有 可 能 结果 的 总 
数 ， 由 整 幅 遥 感 图 像 中 的 端 元 数目 决定 ， 比 如 ， 一 幅 高 光谱 
图 像 的 最 大 端 元 数目 为 3， 那 么 输出 层 的 节点 数 output. 为 
ouput=G+G+G=7. 其中, ANN 的 输出 是 包含 {+1,-1} 的 向 量 ， 
当 ANN 估计 端 元 结果 是 output 中 的 某 一 个 时 ， 其 对 应 的 向 
量 值 为 {t+1}， 同 时 该 像 元 估计 输出 结果 的 其 他 向 量 值 为 {~1}。 
另外 , ANN 的 隐 含 屋 和 输出 层 的 激活 函数 分 别 为 Sigmoid K 
数 和 双 曲 正切 函数 。 

2) 获取 样本 数据 

在 对 ANN 训练 之 前 ， 需 要 准备 样本 数据 ， 包 括 训练 样 
本 和 测试 样本 。 对 于 训练 样本 ， 本 文 首先 随机 设 定 符合 条 件 
的 输出 向 量 ， 然 后 根据 每 一 个 像 元 的 端 元 输出 情况 ， 确 定 参 


LH wr 


x 


出 Xt 


好 的 群 智能 优化 算法 代替 传统 的 梯度 算法 对 目标 函数 进行 求 
解 ， 最 终 实现 高 光谱 图 像 的 解 混 。 


与 混合 的 相关 端 元 ， 再 由 Dirichlet 分 布 [3 产生 丰 度 和 矩阵， 分 
别 基 于 不 同 的 混合 模型 ， 添 加 白 噪 声 进 行 混 合 ， 最 终 获 得 网 
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络 输入 ， 即 训练 样本 ， 同 样 的 方法 也 可 以 获得 测试 样本 。 算法 的 示意 图 。 


3) 训练 网 络 


原 有 经 验 风险 的 基础 上 添 力 


其 定义 为 


在 训练 阶段 ， 本 文采 月 
训练 ， 同 时 为 了 防 J 


昌 反 向 传播 算法 PJ 对 网 络 进行 学 习 
| 练 过 程 中 出 现 过 拟 合 的 情况 ， 以 
及 提高 网 络 的 泛 化 能 力 ， 本 文 引 入 风险 最 小 化 的 策略 ， 即 在 
E 则 化 项 作为 网 络 的 损失 函数 ， 


IMSE = acce -y)-ü- DEX wy) (7) 
m na 


实际 的 网 络 输出 ; 


其 中 : NJ ANN frd 


层 的 节点 数 ; y; 和 si 分 别 为 目标 输出 和 


其 控制 着 向 最 优 误差 寻找 


个 数 。 


3.2 改进 的 差分 搜索 算法 
差分 搜索 算法 是 土耳其 学 者 Civicioglu 根据 生物 体 在 迁 


徙 过 程 中 寻找 有 利生 存 


KER CEP BIO, 提出 的 一 种 群 智能 优 


化 算法 。 文 献 [9] 比 较 了 
题 下 的 优化 性 能 ， 并 证 


DSA 与 其 他 优化 算法 在 相同 求解 问 
其 有 效 性 。 但 是 原始 算法 在 种 群 


的 方向 ， 可 以 加 快 寻 优 的 速度 ; n 为 整个 网 络 所 有 的 参数 的 


搜索 的 过 程 中 ， 由 了 


的 增加 做 相应 的 调整 ， 以 及 当前 个 体 和 寻 优 个 体 之 间 没 有 进 
行 实时 的 参数 控制 ， 导 致 算法 进化 后 期 收敛 速度 下 降 、 局 部 
收敛 精度 降低 ,为 了 平衡 DSA 的 全 
力 ， 加 快 收敛 速度 ， 本 文 在 9 
适应 平衡 矩阵 k， 使 得 种 群 在 进化 


于 搜索 方向 和 搜索 步 长 没有 根据 迭代 次 数 


搜索 能 力 和 局 部 搜索 能 
羊 的 迁徙 方程 中 添加 了 自 


期 具有 较 好 的 全 局 收敛 


的 迁徙 方程 为 


能 力 ， 同 时 在 进化 后 


期 也 有 较 高 的 局 部 收敛 精度 ， 其 改进 后 


StopoverSite = k x Pos + R x map x (donor — Pos) (8) 


其 中 :大 是 一 个 自 适 应 平衡 矩阵 ,=w*m， 其 中 ，w=71-007*25 ， 


epk 为 种 群 当前 迭代 次 数 ，epoch IPREN BEKARE m 
是 一 个 服从 均值 为 1 方差 为 0.3 的 均匀 分 布 产生 的 一 个 大 小 


为 NXD 的 随机 矩阵。 


改进 后 的 搜索 方程 ，w ABS 


适应 调整 。 在 循环 的 初 
力 


从 次 数 的 增加 而 进行 自 


Hj. w 的 值 较 大 ， 算 法 的 全 局 搜索 能 
会 被 放大 ， 随 着 迭代 次 数 的 增加 ，w 的 值 会 慢 慢 变 小 ， 算 


法 的 局 部 搜索 能 力 会 被 变 大 。 本 文采 / 


改进 后 的 DSA 作为 目 


标 函 数 的 优化 方法 ， 通 过 将 丰 度 值 和 待 求 参数 映射 为 种 群 个 


体 的 位 置 变量 ， 可 以 在 少量 的 种 群 个 体 以 及 较 少 的 迭代 次 数 


的 情况 下 有 效 的 实现 解 混 。 


3.3 目标 函数 的 确 


针对 基于 模型 的 钥 


测 值 y EE IR TOES 
接近 ， 于 是 ， 可 确定 解 混 算法 的 
J(à,j) 2 min y- $ IÈ (9) 
个 像素 点 的 丰 度 值 和 模型 中 


3.4 解 混 算 法 步骤 


函数 ， 依 据 最 优化 理 


需要 确定 一 个 合适 的 目标 


谱 图 像 解 混 的 实际 问题 转换 为 


对 目标 函数 的 求解 问 
FREA. KLETE 
KEE a S. H e, Ra 
计 正 确 的 话 ， 那 么 高 光 训 


| 端 元 的 丰 度 值 以 及 模型 的 
EHR ÆRA (reconstruction error, RE) 
度 和 模型 的 未 知 参数 估 
的 图 像 中 单 像素 点 的 观 


tT 


构 的 该 像素 点 高 光谱 图 像 和 撩 量 3 将 非常 


其 中 : 4 和 7 分 别 为 待 估计 的 六 
的 待 求 参数 (LMM 除外 )。 


在 利用 ANN 获得 每 个 混合 像 元 的 端 元 数目 和 类 别 以 后 ， 
后 基于 像 元 混合 模型 确定 各 个 混合 像 元 解 混 的 目标 函数 ， 
行 优化 求解 , 求 得 丰 度 向 


之 后 基 
最 后 采用 改进 的 DSA 对 
量 ， 完 成 对 遥感 图 像 的 解 混 。 


展示 了 本 文 所 提 整 个 解 混 


解 混 算 法 的 具体 实现 过 程 如 下 : 

a) 准 备 训练 样本 。 首 先 利 用 HySime 算法 号 ] 确 定 整 幅 遥 
感 图 像 的 最 大 端 元 数目 ， 并 用 VCA 算法 29 提 取 端 元 光谱 
线 , 由 此 设 定期 望 网 络 输出 的 每 个 像 元 所 有 的 端 元 估计 情况 ， 
同时 采用 Dirichlet 分 布 产生 丰 度 和 矩阵， 根据 网 络 输出 结果 以 
及 像 元 混合 模型 添加 一 定 信 噪 比 的 高 斯 白 噪声 合成 训练 样本 。 

b) 利 用 同样 方法 合成 测试 样本 ， 设 置 并 构建 ANN 网 络 ， 
并 对 网 络 进行 训练 测试 , 最 后 保存 符合 条 件 ( 估 计 准 确 率 95% 
以 上 ) 的 网 络 。 

c) 输 入 高 光谱 遥感 图 像 ， 利 用 训练 好 的 网 络 对 遥感 图 像 
各 个 像 元 的 端 元 数目 和 类 别 进 行 估 计 ， 获 取 单 个 像素 点 的 估 
计 结 果 。 

d) 基 于 式 (9) 确定 各 个 混合 像 元 解 混 的 目标 函数 。 

e) 将 目标 函数 中 的 丰 度 向 量 以 及 待 求 参数 映射 到 改进 的 
DSA 中 的 位 置 向 量 ， 同 时 设 定 DSA WEARER EA 
代 次 数 、 边 界限 制 等 参数 。 

有 初始 化 种 群 位 置 ， 得 到 初始 丰 度 估计 ， 根 据 像 元 混合 
模型 得 到 重 构 的 观测 数据 ， 并 求 得 初始 适应 度 值 。 

2) 将 种 群 按 照 式 (8) 进行 搜索 更 新 ， 并 进行 边界 限制 。 

h) 根 据 DSA 的 贪 焚 选择 机 制 ， 不 断 寻 优 ， 得 到 整个 种 群 
的 全 局 最 优 值 。 

i) 当 达到 迭代 次 数 以 后 ， 输 出 当前 最 优 值 ， 完 成 单个 像 
素 点 的 解 温 ， 否 则 继续 优化 。 


hon 


Endmember estimation based 
on pixels mixing models 


Hyperspectral images Abundance maps 
AINT 2 I£É 
图 2 解 混 算法 示意 图 


Fig.2 Schematic diagram of unmixing algorithm 
4 ”实验 结果 与 分 析 


为 了 验证 本 文 所 提出 的 基于 ANN 端 元 估计 的 高 光谱 图 
像 解 混 算法 的 性 能 ， 接 下 来 将 分 别 使 用 仿真 数据 和 真实 的 遥 
感 数 据 进行 实验 。 本 文采 用 的 解 混 性 能 评价 指标 有 均 方 根 误 
25 27] (root mean square error, RMSE ) 3E H x 26 
(Reconstruction Error, RE), 2615 fü 4) f 05] (spectral angle 
mapper, SAM)。 其 中 ， 重 构 误 差 和 光谱 角 分 布 是 从 算法 重 构 
出 的 观测 数据 和 实际 高 光谱 数据 之 间 的 误差 来 评价 算法 的 性 
能 ， 而 均 方 根 误差 用 来 评价 真实 丰 度 和 算法 估计 的 丰 度 之 间 
的 相似 性 。 三 种 评价 指标 的 定义 如 下 : 
1) RMSE 


ml 


> (10) 


中 :ay 和 si 分 别 为 高 光谱 图 像 第 i 个 端 元 在 第 j 个 像素 点 的 
真实 丰 度 和 估计 丰 度 ; P 入 分别 为 高 光谱 图 像 的 端 元 数目 
和 像素 数目 。 

2) RE 


1 [1 
RMSE = 一 》 | 一 > (o -s,) 
id XN ja 


AL 


IA fig 
RE= 一 》,| 一 > -7 
Lex (=) 


其 中 江 为 高 光谱 图 像 的 波段 数目 ，” 和 > 分 别 为 高 光谱 图 像 
的 第 i 个 像素 点 的 真实 观测 数据 和 重 构 数据 。 
3) SAM 
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S - LY 0b». So) (2) 表 1 标 场景 1 下 各 算法 的 n iod 
p Table 1 Comparison of RMSE values of each algorithm in target scene 1 
"PP — m jn BE SNR-40dB SNR-30dB SNR-20dB SNR=10dB 
s MIN LMM 算法 0.3874 0.3945 0.3838 0.3598 
41 仿真 实验 本 文 算法 0.0042 0.0130 0.0141 0.0434 
仿真 实验 的 高 光谱 数据 来 自 美 国 地 质 勘测 局 (USGS) GBM 算法 0.0221 0.0242 0.0374 0.0681 
提供 的 端 元 光谱 库 ， 库 中 的 光谱 曲线 数据 为 224 个 波段 。 为 本 文 算法 0.0035 0.0095 0.0192 0.0181 
了 使 实验 更 具有 说 服 力 ， 本 文选 择 两 个 目标 场景 的 图 像 进行 PPNMM 算法 0.0232 0.0376 0.0473 0.0842 
解 混 : a) 从 光谱 库 随 机 选取 了 Actinolite, Almandine, Brucite 本 文 算法 0.0019 0.0123 0.0182 0.0559 
这 三 种 端 元 分 别 基 于 LMM、GBM、PPNMM、MLM 4 种 混 MLM 算法 0.0182 0.0276 0.0477 0.0802 
合 模 型 来 合成 50X50 大 小 的 高 光谱 图 像 , 并 在 每 幅 图 像 混合 本 文 算法 0.0045 0.0112 0.0282 0.0359 
的 过 程 中 依次 添加 信 噪 比 (SNR) 为 40dB. 30dB. 20dB. X2 目标 场景 1 下 各 算法 的 RE 值 比较 
10dB 的 白 噪声 ， 其 中 ， 混 合 过 程 中 的 丰 度 矩阵 由 Dirichlet Table2 Comparison of RE values of each algorithm in target scene 1 
分 布 产 生 ， 且 满足 ANC 和 ASC 约束 ; bo 从 光谱 库 随 机 选 RE SNR=40dB SNR=30dB SNR=20dB SNR=10dB 
取 Almandine，Alunite,，Olivine，Andradite，Gaylussite. 这 五 种 LMM 算法 0.0071 0.0223 0.0709 0.2226 
端 元 光谱 按照 第 一 种 场景 的 混合 方式 进行 合成 。 两 种 场景 的 本 文 算法 0.0060 0.0226 0.0632 0.2221 
端 元 光谱 曲线 图 如 图 3、4 所 示 。 GBM 算法 0.0109 0.0249 0.0756 0.2369 
id 本 文 算法 0.0073 0.0178 0.0557 0.0639 
Actinolite NMNHR16485 PPNMM 算法 0.0171 0.0264 0.0725 0.2237 
Almandine WS475 
1.2 Brucite HS247.3B 本 文 算法 0.0054 0.0173 0.0533 0.1734 
MLM 算法 0.0167 0.0254 0.0525 0.1023 
1 本 文 算 法 0.0080 0.0121 0.0330 0.0821 
R3 目标 场景 1 下 各 算法 的 的 SAM 值 比较 
$ Table 3 Comparison of SAM values of each algorithm in target scene 1 
$ 0.6 E SAM SNR-40dB SNR-30dB SNR-20dB SNR=10dB 
= es LMM 算法 0.0100 0.0320 0.1002 0.3074 
9.4 p. 本 文 算法 0.0091 0.0316 0.0995 0.3063 
"M GBM 算法 0.0113 0.0323 0.1000 0.3069 
本 文 算 法 0.0061 0.0184 0.0006 0.0008 
0 PPNMM 算法 0.0254 0.0399 0.1042 0.3236 
0 50 100 150 200 250 本 文 算法 0.0059 0.0192 0.0594 0.1867 
Bands MLM 算法 0.0244 0.0300 0.0712 0.1896 
图 3 目标 场景 1 的 端 元 光谱 曲线 图 本 文 算法 0.0062 0.0100 0.0470 0.0611 
Fig.3 Endmember spectral of target scene 1 表 4 标 场景 2 下 各 算法 的 RMSE 值 比较 


Table 4 Comparison of RMSE values of each algorithm in target scene 2 
Almandine WS475 


Alunite HS295.3B RMSE SNR=40dB SNR=30dB SNR=20dB SNR=10dB 


om ho tm LMM 算 法 —— 04109 0.4447 0.4459 0.4447 

Gaylussite NMNH102876-2 本 文 算法 0.0263 0.0379 0.0544 0.0889 

2 GBM 算法 0.0541 0.0518 0.0542 0.0783 
E 本 文 算法 0.0122 0.0243 0.0369 0.0651 
E PPNMM 算法 ^ 0.3647 0.3654 0.3987 0.4077 
c 本 文 算法 0.0232 0.0339 0.0500 0.0885 
MLM 算法 0.0789 0.0812 0.1032 0.1233 

本 文 算法 0.0521 0.0623 0.0812 0.0911 


表 5 目标 场景 2 下 各 算法 的 RE 值 比较 


Table 5 Comparison of RE values of each algorithm in target scene 2 


0 50 100 150 200 250 RE SNR-40dB SNR-30dB SNR-20dB SNR=10dB 

Bands LMM 算法 0.0048 0.0154 0.0486 0.1543 

图 4 目标 场景 2 的 端 元 光谱 曲线 图 本 文 算法 0.0081 0.0116 0.0330 0.1387 

Fig.4 Endmember spectral of target scene 2 GBM 算法 0.0165 0.0223 0.0532 0.1621 

为 J 验证 基于 ANN 端 元 估计 算法 的 有 效 性 和 抗 噪 , 本 文 算法 0.0096 0.0328 0.0473 0.0240 

本 文 分 别 基 于 LMM, GBM, PPNMM, MIM 在 同样 的 条 件 PPNMM 算法 0.0084 0.0165 0.0498 0.1535 

了 两 种 目标 场景 下 使 用 原始 算法 和 本 文 算法 在 不 同 信 噪 比 条 MIME 00091. 0003 — 0031 — 00D1 

件 下 的 RMSE, RE, SAM 值 。 为 了 避免 随机 性 ， 所 有 实验 本 文 算法 0.0078 0.0091 0.0123 0.0334 
结果 均 为 独立 进行 10 次 以 后 所 取 的 平均 值 。 
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表 6 目标 场景 2 下 各 算法 的 SAM 值 比较 


Table6 Comparison of SAM values of each algorithm in target scene 2 


SAM SNR=40dB SNR=30dB SNR=20dB SNR=10dB 
LMM 算法 0.0321 0.0321 0.1013 0.3110 
本 文 算法 0.0210 0.0210 0.0845 0.2150 
GBM 算法 0.0376 0.0376 0.1032 0.3113 
本 文 算法 0.0210 0.0210 0.0847 0.0317 
PPNMM 算法 0.0348 0.0348 0.1046 0.3143 
本 文 算法 0.0227 0.0227 0.0769 0.2492 
MLM 算法 0.0224 0.0224 0.0490 0.0812 图 5 Samson 三 维 立 体彩 色 图 
本 文 算法 0.0107 0.0107 0.0212 0.0512 Fig.5 Samson 3D color map 
实验 表格 可 以 发 现 : 


a) 无 论 是 目标 场景 1 还 是 目标 场景 2， 基 于 不 同 的 像 元 
混合 模型 ， 在 相同 的 条 件 下 ， 利 用 ANN 选择 每 个 像素 的 相 LMM 算法 
关 端 元 进行 解 混 的 精度 比 解 混 过 程 中 保持 端 元 数目 不 变 的 算 
法 精度 更 高 。 可 以 验证 ANN 端 元 估计 的 方法 是 有 效 的 ， 而 
结果 是 相对 稳定 的 。 

b) 将 表 中 不 同 SNR 下 的 RMSE、RE、SAM 值 对 比 可 以 
AOL, 所 有 解 混 算 法 的 性 能 都 会 随 着 SNR 的 减 小 而 逐渐 变 差 ， GBM 算法 
而 且 原 有 模型 算法 的 实验 结果 比 本 文 算法 的 结果 下 降 的 更 加 
， 可 以 发 现 ， 基 于 ANN 端 元 估计 解 混 算法 的 抗 噪声 性 
中 鲁 棒 性 均 优 于 传统 的 解 混 算法 。 此 外 ， 在 相同 的 实验 条 
EF, GBM 作为 一 个 广义 通用 的 解 混 模 型 具有 最 佳 的 解 混 
性 能 。 PPNMM 算法 
c) 将 两 个 目标 场景 的 结果 对 比 发 现 , 当 端 元 数目 增加 时 ， 
于 估计 输出 结果 的 增加 以 及 ANN 结构 自身 的 局 限 性 ， 本 
文 解 混 算 法 的 性 能 有 所 下 降 ， 但 是 总 体 上 还 是 高 于 现 有 解 混 
算法 的 效果 。 因 此 ， 本 文 提出 的 利用 ANN 估计 端 元 的 解 混 
m 策略 是 有 效 的 ， 同 时 本 文 所 用 ANN 结构 简单 ， 在 获取 每 个 MLM 算法 
AN] 像 元 的 端 元 数目 和 类 别 以 后 ， 可 以 实时 减少 模型 参数 ， 减 小 
re 计算 量 ， 和 传统 的 解 混 算法 相 比 ， 具 有 更 大 的 优越 性 。 
4.2 真实 遥感 数据 实验 
为 了 进一步 验证 本 文 所 提 的 基于 ANN 端 元 估计 的 高 光 
(VY 谱 图 像 解 混 算法 的 性 能 ， 选 取 由 美国 Florida Environmental 本 文 算法 
pm Research Institute 拍摄 的 Samson 数据 作为 高 光谱 解 混 的 测试 


b x 
GG 月 
Lili 
[und 


iS 


对 象 ， 该 数据 中 一 共 记录 了 156 个 波段 ， 其 波段 范围 为 E" 
401-889 nm,， 光 谱 分 辩 率 为 3.13 nm， 由 于 原始 图 像 共 有 952 (a) 岩 石 (b) 水 (OPR 
X952 个 像素 点 ， 图 像 较 大 ， 为 了 降低 计算 复杂 度 ， 通 常 选 图 6 各 种 算法 的 解 混 丰 度 图 

取 如 图 5 所 示 的 95x95 子 区 域 进行 实验 ,该 图 像 区 域 主 要 包 Fig.6 Unmixing abundance maps of various algorithms 
含 树木 、 岩 石和 水 三 种 地 物 。 表 7 不 同 解 混 算 法 的 结果 比较 

在 本 节 中 ， 为 了 更 加 直观 地 验证 本 文 算 法 的 性 能 ， 同 时 Table 7 Comparison of results of different unmixing algorithms 

鉴于 GBM 具有 广泛 的 适应 性 和 较 大 的 灵活 性 ， 既 能 详细 地 比较 算法 RE SAM 
雪 述 非 线性 混合 的 地 物 场景 ， 也 能 涵盖 线性 混合 场景 ， 本 文 LMM 算法 0.0990 0.1723 
选 定 GBM 作为 所 提 解 混 算 法 的 像 元 混合 模型 ， 并 用 实际 的 GBM 算法 0.0754 0.1547 
高 光谱 遥感 图 像 数据 进行 验证 。 由 于 实际 场景 的 先 验 信息 缺 PPNMM 算法 0.0670 0.1561 
失 ， 本 文 首先 使 用 VCA 算法 获取 整 幅 遥感 图 像 的 所 有 端 元 MLM 算法 0.0660 0.1449 
的 光谱 曲线 ， 之 后 基于 GBM 模型 并 添加 白 噪声 随机 合成 训 本 文 算法 0.0321 0.0702 
练 样本 来 训练 网 络 ， 并 用 符合 条 件 的 网 络 估计 实际 遥感 图 像 5 ”结束 语 
各 像素 点 的 端 元 数目 和 类 别 ， 最 终 采 用 改进 的 差分 搜索 算法 SPESE 


对 目标 函数 进行 优化 求解 ， 完 成 解 混 。 表 7 给 出 了 本 文 算法 本 文 提出 一 种 基于 ANN 端 元 估计 的 高 光谱 图 像 解 混 算 
与 LMM SEU, GBM 算法 [031]、PPNMM 算法 09、MLM 算 法 ， 首 先 利 用 ANN 对 遥感 图 像 各 个 像 元 的 端 元 数目 和 类 别 
ikUSIB] RE, SAM 值 ， 图 6 给 出 了 几 种 算法 的 实际 解 混 丰 度 ”进行 估计 ， 之 后 根据 估计 结果 确定 目标 函数 将 解 混 问题 转换 
图 。 为 目标 函数 优化 问题 ， 同 时 利用 改进 的 差分 搜索 算法 对 目标 
根据 实验 结果 以 及 对 比 实 际 的 解 混 丰 度 图 可 以 发 现 ， 本 ”函数 进行 优化 求解 ， 最 终 实现 图 像 解 混 。 仿 真实 验 和 真实 的 
文 基于 ANN 端 元 估计 的 高 光谱 图 像 解 混 算 法 和 现 有 主流 的 遥感 数据 实验 表明 , 本文 所 提 解 混 算 法 不 仅 提高 了 解 混 精 度 ， 
解 混 算法 相 比 ， 具 有 更 高 的 解 混 精 度 和 更 优 的 解 混 效果 。 而 且 克 服 了 梯度 求解 算法 的 局 限 性 ， 最 终 改善 了 解 混 算法 的 
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效果 。 同 时 ， 本 文选 取 的 ANN 结构 简单 ， 对 于 地 物种 类 数 
目 较 多 的 场景 , 解 混 算法 的 性 能 有 所 下 降 ,在 今后 的 研究 中 ， 
可 采用 更 加 复杂 的 网 络 模 型 比如 卷 积 神经 网 络 
(Convolutional Neural Network, CNN) BI、 深度 置信 网 络 
(Deep Belief Network, DBN) B0 等 ， 进 一 步 提 高 解 混 算法 的 


性 能 。 
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